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1. INTRODUCTION 

The fracture strength of composite materials in 
general and ceramics in particular depends to a consid- 
erable extent on the size, shape, orientation, and 
distribution of flaws or imperfections in the material. 

In studying the "bulk strength" and other bulk properties 
of the composite it is generally assumed that these flaws 
are randomly distributed throughout the material and 
hence the medium is statistically homogeneous. There- 
fore, in this type of studies the very nature of the 
phenomenon requires that some kind of a statistical 
strength theory be used as a guide in the investigations. 

In a second type of approach to fracture studies, 
one is basically interested in the more detailed treat- 
ment of the fracture initiation phenomenon in the local- 
ized regions. In such an approach it is generally 
conjectured that the macroscopic fracture of the medium 
will initiate at and will propagate from a "dominant 
flaw". This is usually the flaw which has, from the 
viewpoint of the fracture resistance of the material, 
the worst possible combination of size, shape, orienta- 
tion, and location in the structural part which has a 
given overall geometry and is subjected to a given 
system of external loads and environmental conditions. 
The dominant flaw may be caused by residual stresses or 
by some other type of loading before the part is put 
into use, it may be a manufacturing flaw (such as en- 
trapped gas, weak impurities, imperfect bonding, and the 
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geometric singularities arising from the particular shape 
of the constituent materials) , or it may result from the 
growth of a "microflaw" due to the cyclic nature of the 
operating stresses. In some cases it may be possible to 
detect such flaws by using nondestructive testing tech- 
niques. However, more often, in studies relating to 
structural integrity and reliability one simply assumes 
their existence. 

In ceramics one may assume that the composite medium 
consists of the matrix, the inclusions with elastic prop- 
erties different than that of the matrix, and the pores. 
Generally, from the viewpoint of the fracture strength, 
the worst imperfections are the cracks in the matrix 
and/or in the inclusions. Thus, in the most elementary 
sense the main objective of fracture studies in ceramics 
is to develop a technique for the solution of the follow- 
ing mechanics problem: To find the magnitude of the 

external loads which can be carried safely by a struc- 
tural part with given matrix-dominant crack-inclusion- 
pore geometry under given environmental conditions . It 
is clear that the problem has three different aspects 
which may be studied more or less independently: 

(a) Determination of the "characteristic strength 
parameter" of the material. This is the experimental 
work aimed at determining the inherent resistance of the 
constituent materials to fracture corresponding to the 
conjectured mode of failure, such as "brittle" or "quasi- 
brittle" (i.e., low energy-type) fracture, "ductile" 

(i.e., high energy-type)’ fracture, shear or tensile 
cleavage, and debonding or delamination. Example: frac- 
ture toughness, G^. 

(b) The solution of the mechanics problem for the 
given geometry and the external loads in order to deter- 
mine the "critical load factor" representing the severity 
of the applied loads and corresponding to the conjectured 
mode of fracture. Example: strain energy release rate, 

G I* 

(c) Development of an appropriate fracture theory 
and a related quantitative "fracture criterion" . This 
is usually nothing but the direct comparison of the 
"characteristic strength parameter" of the material with 
the "critical load factor". Example: Gj<Gxc* It is 
clear that in fracture studies the development of the 
fracture theory and selection of a proper fracture 
criterion must precede the experimental determination of 
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the characteristic strength parameters and the theoret- 
ical evaluation of the corresponding critical load 
factors. 

In this paper after a brief review of some of the 
current fracture theories, a group of mechanics problems 
of practical interest relating to the elastic interaction 
between cracks and inclusions will be identified and a 
summary of the results will be presented. 


2 . FRACTURE THEORIES 

Due to the lack of full physical understanding of 
the fracture phenomenon in solids and the lack of suffi- 
ciently powerful mathematical tools, at the present time 
there is no consistent theory dealing with all the rele- 
vant aspects of the phenomenon. Partly for this reason 
and partly because of the divergence in the background 
and interest of 'the researchers studying the field, 
generally the existing theories treat the subject from 
only one of the three points of view, namely (a) atomic 
(or microscopic) , (b) microstructural , and (c) continuum. 

From the viewpoint of structural applications the theo- 
ries which provide the quantitative tools, to deal with 
the problem are mostly the continuum theories which are 
based on the notions and principles of the classical 
thermodynamics and the continuum solid mechanics. Hence, 
the discussion given in this section will include only 
these theories. 


(a) Energy Balance Theories 

The fracture theories which are based on the energy 
balance follow directly from the thermodynamic principle 
which states that in a given continuum 

U* = t’ + v’ + D* , (1) 

where U is the work done by the external loads, T is the 
kinetic energy, V is the (recoverable) potential energy, 
D is the total dissipation, and the dot stands for the 
time rate of change [ 1 ] . (1) is valid for the whole as 
well as a toroidal portion of the body enclosing the 
leading edge or the periphery of the crack. At the on- 
set of (quasi-static) crack propagation in ideally 
brittle materials (1) is easily shown to reduce to the 
Griffith criterion, namely 
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G IC 2y s ' 


( 2 ) 


= 


where Gj is the "strain energy release rate" or the 
"crack extension force" per unit length of the crack 
front, and y s is the specific surface tension energy. 

In the extension of the theory to quasi-brittle type 
fracture [2,3] Gjc is interpreted as the "fracture 
toughness" of the material, or Ys i s replaced by Yp, the 
"fracture energy" (i.e., the amount of energy required 
to create a unit area of fracture surface) . In either 
case, noting that Gj = , the criterion may 

also be stated in terms of the "stress intensity factor" 
as follows: 


K 


I 


K IC 


(3) 


where Kj is the stress intensity factor (for example 
K I = ^l/i" = a 0 /ifa for an infinite solid with a through 
crack of length 2a and subjected to a uniform stress 
perpendicular to the plane of the crack) , and K IC is 
assumed to be a material constant. 


Even' though the underlying principle of the theory 
is quite general, quantitatively it gives good results 
only for "plane strain" type of configurations where the 
characteristic size of the "dissipation zone" surrounding 
the leading edge of the crack is small in comparison with 
the characteristic crack size. Thus, the subscript I 
in (2) and (3) is usually meant to stand for "plane 
strain". In the general case writing [2,3] G=G C , one 
may interpret the load factor G as the "energy available" 
(through external addition or internal release) for 
fracture, and the characteristic strength parameter Gc 
as the "energy required for the creation of unit fracture 
area". In practice the difficulty lies in evaluating G 
with sufficient accuracy and in properly isolating and 
measuring G^. The present chaotic state of characteriza- 
tion of the so-called "plane stress fracture" is a good 
example. Even though developed only for nonlinear elas- 
tic (hence, nondissipative) solids, recent attempts of 
computing G through path-independent integrals (known as 
the J-integral) have met with partial success [4-6]. 
However, the problem of a reliable characterization for 
G q still remains very elusive.* It should be noted that 

*The recent investigations on the development of the 
"crack extension resistance curve" for (or the so-called 
Kj^-characterization of) the fracture of sheet materials 
appear to be promising. However, in this case too the 
isolation of a geometry-independent characteristic 
strength parameter has not yet been successful. 
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regardless of the generality of its underlying principle, 
as it is understood and applied in practice, a tacit 
requirement of the energy balance theory is that the 
stress and deformation fields remain "self-similar" as 
the crack propagates, which in a narrow sense means that 
the crack remains in its ov/n plane as it propagates, the 
crack front advancing in an autonomous fashion. This is 
the case only if {a) the crack remains imbedded in the 
same homogeneous medium as it propagates, i.e., it does 
not intersect any phase boundary or bi-material inter- 
face, and (b) the crack is and remains as an interface 
crack. For example if ki and k 2 are the in-plane tensile 
and shear components of the stress intensity factor* at 
the tip of a crack lying along the interface of two 
elastic materials with constants y]_,Kp and \±2' K 2 (i-i : 

shear modulus, Ki = 3-4 Vi for plane strain, Kj_ = (3-Vi)/ 
(1+vi) for plane stress, Vi: Poisson's ratio), the 
(elastic) strain energy release rate (per unit crack 
extension/per unit crack front) is given by [7] 


G 


which, 



■n (y l + <ih 2 ) ^2 + ic 2 y l ) „ 2 2, 

4 y 1 y 2 [ (1 +k 1 )u 2 + (l+K^yJ^l K 2 > ’ 


for homogeneous materials reduces to 


(4) 


... 3_U _ TT (1 + K) 
3a 8y 


(k 1 2 +k 2 2 ) . 


(5) 


In (4) and (5) , again the assumption is that the crack 
remains in its own plane as it propagates . In applying 
(4) the related characteristic strength parameter, (Gc )^2 
will be the fracture energy of the particular adhesion. 


(b) C.O.D. -Characterization 

Contending that fracture is basically creation of 
new surfaces in the solid resulting from rupture- 
separation of the material around the leading edge of 
the crack, particularly in the presence of large scale 
inelastic deformations, some investigators have concen- 
trated their attention at the crack tip region and pro- 
posed the argument that a true measure of the severity 


*In this paper the definition of the stress intensity 
factors k-j_,k 2 ,k 3 does not include /if and is based on 
the asymptotic relations of the form crgg(r,0) - k^/Z^F , 
a r g(r,0) - k 2 //2r , crg z (r,0) = k 3 // 2 r for small r, r 
and 6 being the polar coordinates at the crack tip. 
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of the external loads may be the deformation state at 
the crack tip characterized by the amount of stretch (<5 
\ or C.O.D., the proposed load factor) suffered by a thin 

layer of material in the plane and ahead of the crack. 

In turn the fracture resistance of the material may be 
characterized by the value of this stretch at rupture 
(5 q/ the proposed strength parameter) which is conjec- 
> tured to be (with certain qualifications)' a material 

{ constant [8-11]. The criterion (which has been standard- 

j ized in the United Kingdom and has found so far a wider 

i acceptability in Europe than in the U.S.) appears to 

5 work rather successfully for "plane stress" type fracture 

: of thin sheets with relatively large plastic zones around 

l the crack. The application of this criterion too re- 

* quires that the extension of the crack be co-planar with 

I the crack and the crack remain imbedded in the same 

l homogeneous medium. 

i . 

* 

| (c) Critical Stress' Criterion 

j ' 

l In ceramics and other composites, because of the 

| nonhomogeneous nature of the medium, usually the stress 

’ state around the leading edge of the crack is of the 

? mixed-mode type, and one occasionally encounters the 

problem of a crack intersecting a bi-material interface. 

I In cases such as these either the crack may not remain 

I co-planar or the stress and displacement fields around 

\ the crack tip may not remain self-similar as the crack 

f propagates.* Thus, in this type of fracture studies the 

I fracture theories and criteria based on the notions of 

:> energy balance or critical crack opening stretch men- 

< tioned above will be either ineffective or inapplicable, 

j and clearly a more flexible approach is needed. For 

brittle and most quasi-brittle materials the "maximum 
■* cleavage stress" criterion discussed in [13] or some 

! variation of it may provide an acceptable and a very 

simple working tool. The modification of the criterion 
is necessary primarily to be able to deal with problems 
in which instead of the standard r -J_ / 2 singularity the 
stresses have either oscillating and/or arbitrary power 
\ singularity (i.e., r~ a , 0<Re(a)<l), as in the case of 

;! cracks lying on or intersecting the interfaces. The 


*As will be pointed out later in this paper for a crack 
tip terminating at an interface the power of the singu- 
larity is no longer -1/2 and the r, 6 -dependence of the 
stresses is completely different than that of a crack 
tip imbedded in a homogeneous medium (see [12]) . 





argument here is as follows: 

Due to the very high concentration of stresses and 
strains and to highly nonlinear and inelastic nature of 
the phenomenon, even in brittle materials, the response 
of the medium in a small toroidal region enclosing the 
leading edge of the crack cannot be accurately obtained. 
Noting that the initial radius of curvature of the crack 
surface at the tip region may be as small as the inter- 
atomic distances, even if the necessary mathematics were 
available, there is always a region in which the contin- 
uum theories of solids will be inapplicable. Let o D be 
the radius of the cross section of this region (measured 
from the edge of the crack) . 6p, called the "character- 

istic process zone size", may depend on the microstruc- 
ture and continuum properties of the material as well as 
on the environmental conditions. Even though what goes 
on in this region (as far as the material response is 
concerned) is very important and may affect the response 
outside the region, it cannot quantitatively be incor- 
porated into a continuum model. Hence any continuum 
theory or criterion will have to be based on the quanti- 
ties which can be calculated and measured outside this 
region. Thus, we are essentially forced to ignore the 
toroidal region (r < 6 p , -rr < 6 < ir) , where r and 6 are the 
polar coordinates at the crack tip in the plane perpen- 
dicular to the crack front. The "critical stress" 
criterion may then be stated as 

the fracture propagation will take place in 
the direction 0 = 0c for which the cleavage 
stress crQ 0 (6p,0 c ) is maximum and when 

°ee ({ p ,9 c ) - °c ■ (6) 

where c c is constant for a given material 
under given environmental conditions . 


(6) is valid if the "weak link" is a particular 
homogeneous medium surrounding the crack tip. However 
if the singular point is at a bi-material interface and 
if the weak link is the adhesion along the interface, 
then (6) may be modified as 


, 2,2 . 1/2 

a 08 + a r0 

a r8 + fa 88 


a 88 >0 / 


' a 00 <0 ' 


(r= <« p ) 12 , e = e 12 ) 


where (a^)i2 represents the adhesive strength of and 
(<Sp)i 2 is the characteristic process zone size for the 
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given bi-^material interface, and 0]_2 is the radial direc- 
tion along the interface. In (7) for Of . q > 0 the right- 
hand side represents the magnitude of the (tensile) 
stress vector acting on the surface of adhesion. For 
000 < 0 the mode of adhesive rupture has to be shear, and 
the shear stress a r g will be reduced somewhat by "fric- 
tion". The term f 00 Q is therefore added to a r Q to 
account for this effect. Finally, in order to be able 
to apply the criterion in practical problems, it will 
be assumed that due to relative "smallness" of 6p 

the effect of the phenomena taking place within 
the process zone r < <$p on the stresses outside 
this region is negligible. 


Thus, for example, on a crack front imbedded in a 
homogeneous brittle medium, if kj_ and k 2 are the tensile 
and shear components of the stress intensity factor in 
a plane perpendicular to the leading edge of the crack, 
we have [13] 


and 

or 


°ee ( V 0) “ 


/26" 


6 .. 2 6 3 

cos^ (k^cos 2 ~ 2 


. . 6 , 
k^in^) 


( 8 ) 


8 C is obtained from (3o 00 / 30) = 0, (3 2 a e9 /36 2 ) < 0 

k 2 (1 - 3cos6 c ) - k^^sinO^, = 0 , 

3k„sin0 - k,cos0_ < 0 . (9.a,-b) 


r 


For a crack terminating' at or intersecting the interface 
the stresses around the singular point at the interface 
are of form (for details see [12] and [14]) 


a i j (r , 0 ) 


a — [k, f, . . (0) 

r 0t L 1 111 


k 2 f 2ij (0 ^ 


(i,j = r,0; 0 < Re (a) < 1 ; -tt<0<tt). (10) 

It is seen that unlike the energy balance and the criti- 
cal C.O.D. criteria, the criterion given by (6) or (7) 
can be applied to fracture problems represented by (8) 
or (10) with equal facility. It should be pointed out 
that for symmetrically loaded (i.e., k 2 = 0 ) homogeneous 
materials, in the absence of large scale inelastic 
deformations around the crack front, all three theories 
described in this section are equivalent and would give 
the same results. 





3. INTERACTION BETWEEN INCLUSION AND CRACKS 

Assuming that we have an acceptable "fracture cri- 
terion" (e.g. , any one of the three criteria discussed 
in the previous section or some other criterion) and 
could perform the idealized experiments to determine the 
related "characteristic strength parameters", the problem 
is then reduced to the evaluation of the "critical load 
factor" for the particular structural part with given 
matrix-crack-inclusion-pore geometry and subjected to a 
given set of external loads and environmental conditions. 
The exact treatment of the three-dimensional problem with 
a regular or a random array of inclusions imbedded into 
a matrix containing a (dominant) internal crack is hope- 
lessly complicated. Therefore, we will first assume that 
the composite medium is basically linearly elastic and 
the effect of the small process or the damage zone around 
the crack on the stress distribution is negligible. 

Hence, in the application of any of the fracture cri- 
teria, to evaluate the related "critical load factor" it 
is sufficient to calculate the stress intensity factors 
ki and k 2 / the power of the singularity a, and the func- 
tions fki j , (k,i,j=l, 2) shown in (8) and (10). The 

remaining assumptions which are made in this paper con- 
cern the idealization of the geometry of the medium in 
order to make the analysis tractable. 

First it will be assumed that the medium is a large 
composite "plate" or "cylinder” which is under either 
in-plane stress state or anti-plane shear. It will then 
be assumed that the elastic matrix contains only sparsely 
distributed inclusions, which are circular. Hence the 
mechanical interaction will be primarily between an iso- 
lated inclusion (or a hole) and line crack (s) located in 
the matrix and/or in the inclusion. As in many problems 
in linear fracture mechanics, the solution of this prob- 
lem too may be obtained through superposition of two 
solutions: that of the uncracked medium under a given 
set of external loads, and that of the cracked medium 
(i.e., the perturbed solution) in which the crack surface 
tractions obtained from the previous solution are the 
only external loads. Clearly, only the second solution 
has the singular stresses. Special cases of this problem 
were considered in [15-17]. The results given in this 
paper are taken from [18-20]. 
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(a) Symmetrically Loaded Radial Cracks 

The problem is described in Figure 1. All the prob- 
lems which are discussed in this paper may be formulated 
in terms of a system of singular integral equations of 
the following form: 


b l f, (t) b l b 2 

/ — - t ~ - dt + / k^ (x, t) f 1 (t)dt + / k 12 (x,t) f 2 (t)dt 

a l a l a 2 

7T (1+K ) 

= — p 1 (x) , (a 1 _<x<b 1 ) / 


b 2 f (t) b l b 2 

/ ~ t - — - dt + / k 21 (x,t) f 1 (t)dt + / k 22 (x,t) f 2 (t) dt 

a 2 a l a 2 

7T(1 + K ) 

= 2ir; p 2 (x) ' (a 2 <x<b 2 ) , (11) 


where pi and p2 are the known crack surface tractions, 
k-Lj(x,t) / (i,j = l,2) are known functions (bounded if the 
cracks are fully imbedded in a homogeneous medium, of 
the generalized Cauchy type if one or both cracks inter- 
sect the interface), and the unknown density functions 
fl and f 2 are the derivatives of the crack surface dis- 
placements (for complete details and the solution see 
[18-20]). 


Even though the formulation is quite general, the 
numerical results are given for only two typical and 
somewhat extreme cases, namely the case of a pore (i.e., 
P 2 = 0) and a hard inclusion with (j^/ui) = 23, = 0.35, 

V2=0.3, which approximately corresponds to a metallic 
inclusion in an epoxy-type matrix. Table 1 shows the 
stress intensity factors for a crack in the matrix (i.e., 
a 2 = ^2' R<a]_<b]_, Figure 1) under uniaxial stress o Q at 
infinity, kq j , (j=l,2) shown in the table are the 

dimensionless quantities given by 


‘11 


°o /5 r 


‘12 


k 1 (a^) 
a o /5 l 


( 12 ) 


where 2cq is the crack length. Note that in the case 
of a hole as a^-*R, k^ 2 ‘ >co - in the inclusion case the 
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Table 1. The stress intensity factors 
for a symmetrically located radial crack 
(C 2 = 0 , R/c^ = 2, k^ = 1.6, i<2=1-8). 



^2 

= 0 

y 2 2 

3y l 

b/c 1 

k ll 

k 12 

k ll 

k 12 

3.0 

3.2 

1.417 

CO 

2.274 

0.827 

-> 0 
0.467 

3.5 

1.290 

1.722 

0.874 

0.671 

4 

1.188 

1.394 

0.918 

0.321 

5 

1.102 

1.174 

0.957 

0.924 

6 

1.065 

1.099 

0.973 

0.959 

8 

1.033 

1.045 

0.987 

0.982 

CO 

1.0 

1.0 

1.0 

1.0 


stress intensity factor ratios decrease as the crack 
approaches the interface, and for a]_ = R the power of 
the singularity becomes a where a = 1/2 for homogeneous 
material 0 < a < 1/2, for y 2 > ^1' anc ^ 1/2 < a < 1 for ^2 < l 1 ! 
[12]. In the present case {i.e., y 2 =23y;] L , k^ = 1.6, 

<2 =1-8) , a is found to be 


a = 0.33811 , for plane strain, 

a = 0.28901 , for plane stress. 


(13) 


In this case the stress 'intensity factor at a^ = R is 
defined by and found to be 

kf (a^) = lim /2 (a, - x) a o 2 (x, 0) = 1.375 a 0 c., a , 


where the subscripts 1 and 2 refer to the matrix and the 
inclusion, respectively. The functions of 0 giving the 
angular distribution of the stresses around x = aj_ = R 
may be found in [12], For the crack terminating at the 
interface, the effect of the modulus ratio l^/^l on a 
and k(a^) is shown in Table 2. Table 3 shows the effect 
of the modulus ratio and the Poisson's ratios on the 
stress intensity factor in a completely cracked inclu- 
sion (i.e., a 2 = -R, b 2 = R» a-^ =b^_) where 


k(R) = lim /2 (x - R) a cr, (x,0) = k'a R a . (15) 

x-*R ^ 


-n- 



Table 2. The effect of modulus ratio on 
the stress intensity factors for a crack 
terminating at the interface (ai = R, 
b- I _/R = 2 / Kj_= <2 = 1-8, c^= (b 1 -a 1 )/2). 


y 2 

a 

k(b 1 ) 

k (R) 

y l 


a 

a o c l 

0 

0.05 

0.81730 

2.808 

1.615 

1.053 

1/3 

0.62049 

1.229 

0.5836 

1.0 

0.5 

1.000 

1.000 

3.0 

0.40074 

0.8610 

1.299 

10.0 

0.33277 

0.7969 

1.389 

23.0 

0.30959 

0.7796 

1.375 

100 

0.29387 

0.7691 

1.345 

300 

0.28883 

0.7667 

1.348 


Table 3. Stress intensity factor for a completely 
cracked inclusion (a^=b^, a 2 =-R, b 2 =R, k ' = k (R) / (c 0 R a ) ) . 


y 2 

<1 < 2 

ll 

H 

• 

CO 

<^= 2.2, 

k 2 =1.8 

<-^=1.8, 

<2 2.2 

<1 < 2 = 

= 2.2 

| 

a 

k' 

a 

k' 

a 

mm 

a 

mm 

- 

0.2 

0.3662 

0.789 

0.3809 

0.785 

0.3203 

1.046 

0.3385 

1.010 

0.6 

0.4503 

1.014 

0.4703 

0.946 

0.4212 

1.174 ! 

0.4447 

1.068 

1.0 

0.5 

1.0 

0.5199 

0.921 

0.4772 

1.107 

0.5 

1.0 

2.0 

0.5745 

0.884 

0.5919 

0.817 

0.5569 

0.947 

0.5762 

0.861 

5.0 

0.6789 

0.656 

0.6912 

0.619 

0.6638 

0.694 

0.6773 

0.650 


Table 4 shows the result for a crack intersecting 
the interface. Here note that there are three singular 
points and the stress intensity factors ky and k X y at 
the singular point x=R are defined by 

k (R) = lim y a a. (R,y) , k (R) = lim y a a, (R,y), 
y y->0 yy y y-*-0 y 

(16) 

and 2c is the total crack length. Further results and 
examples showing the crack surface displacements may be 
found in [18] . 





























Table 4. The stress intensity factors for a crack 
crossing the interface (k-|_=1.6, k 2 = 1.8, 

M 2 /Mi= 23 .077 , a=0. 27326 , c = (b 1 -a 2 ) /2) . 


k(b x ) 


c 

o 


0.548 

0.513 

0.570 

0.626 

0.672 

0.710 

0.765 

0.811 


0.920 

0.954 

0.757 

0.670 

0.598 

0.547 

0.518 

0.510 

0.525 

0.572 


0.487 

0.425 

0.619 

0.731 


k(a 2 ) 


<r /c 
o 


3.564 

3.701 

3.756 

3.799 

3.838 

3.874 

3.935 

3.996 


5.600 
5.242 
5.003 
4.570 
4.034 
3.481 
2.956 
2.465 
1.981 
1.383 
-*■ 0 

3.377 
2.978 
1.961 
1.401 
-*• 0 



■0.847 

■0.446 

■0.282 

■0.171 

■0.0335 

-0.0113 

0.105 

0.219 


■ 1.102 

■0.954 

■0.753 

-0.547 

■0.363 

-0.209 

-0.0852 

0.0171 

0.108 

0.212 


■0.730 

■0.229 

0.211 

0.412 


0.170 

0.0894 

0.0565 

0.0342 

0.0167 

0.00227 

- 0.0211 

-0.0440 


0.221 

0.191 

0.151 

0.110 

0.0727 

0.0420 

0.0171 

-0.00343 

-0.0216 

-0.0426 


0.146 

0.0459 

■0.0424 

■0.0826 


(b) Plane Problem for Arbitrarily Located Crack 

The results for the general problem of an arbitrar- 
ily located crack in a matrix containing an inclusion 
or a hole are shown in Tables 5-9. Table 5 shows the 
results for the case in which the crack and the inclusion 
are symmetrically located with respect to the y-axis 
(i.e., for b = 0). Here c is the distance of the crack 
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Table 5. Stress intensity factors for a symmet- 
rically located "tangential crack" perpendicular 
to the load (b=0, R = 2a, <-^=1.6, <2=1-8). 



V 2 /Vj 

= 0 

v 2 /v 1 

= 23 

c/a 

k ll =k 12 

k 21 =-k 22 

k ll =k 12 

k 21 =-k 22 

2.1 

0.0626 

-0.170 

1.201 

0.105 

2.5 

0.260 

-0.101 

1.196 

0.171 

3 

0.395 

-0.0644 

1.196 

0.202 

4 

0.574 

-0.0175 

1.180 

0.180 

6 

0.768 

0.0274 

1.110 

0.109 

10 

0.906 

0.0361 

1.046 

0.0567 

CO 

1.0 

0 

1.0 

0 


from the x-axis (for the notation see the inserts in 
Figures 2 and 3) . In Tables 5-9 the normalized stress 
intensity factors are defined by 


i 



k^ (b+a) 

a /a 
o 


I 



k 2 (b+a) 


(17) 


k 



k x (b-a) 

a /S' 
o 



k 2 (b-a) 

a / a 

o 




i.e., the first subscripts 1 and 2 respectively refer 
to the first and second mode, and the second subscripts 
refer to the right and left crack tips. Again in all 
the results given in Tables 5-9 the material constants 
are P2 =0 ' or Kq=l-6, <2 = 1.8, and the 

external load a Q is perpendicular to the plane of the 
crack. From Table 5 it is seen that k^ may be greater 
than 1 for the stiffer inclusion and considerably less 
than 1 for the hole. This and the similar results ob- 
served in the subsequent tables may at first appear to 
be somewhat paradoxical. However if one considers the 
distribution of the stresses in an uncracked matrix 
shown in Figures 2 and 3, the explanation for the trends 
would be apparent. 

Tables 6-9 show some further results for an arbi- 
trarily oriented crack in the neighborhood of an elastic 
inclusion or a hole. Note that in these examples the 
shear component of the stress intensity factor k 2 j , 
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(j=l/2) is not zero. Thus from (9) it can be shown that 
at the crack tip x=b~a (see Figures 2 and 3) the direc- 
tion of probable crack propagation will be towards the • 
interface if y 2 = 0 (or generally if P2 < U]_) , and away 
from the interface if P2 > Pi- The details of the arbi- 
trarily oriented crack problem and further results may 
be found in [19]. 


Table 6 . Stress intensity factors for a crack perpendi- 
cular to the external load (R=2a, b=3a, <■> =1.6, 

< 2 = 1 . 8 ) . 


2 - 0 m 2 - 23y 1 



c/a 

k ll 

k 12 

k 21 

CN 

CM 

k ll 

k 12 

k 21 

k 22 

0 

0.3 

1.607 

-> 00 
4.267 

-> 0 
0.072 

-> 0 

-0.391 

0.784 

-»■ 0 
0.225 

0 

-0.004 

+ 0 
0.072 

0.5 

1.544 

3.070 

0.113 

-0.412 

0.792 

0.341 

-0.006 

0.101 

1.0 

1.430 

1.969 

0.177 

-0.316 

0.817 

0.613 

-0.005 

0.057 

1.5 

1.357 

1.552 

0.184 

-0.241 

0.839 

0.763 

0.008 

-0.007 

2 

1.299 

1.337 

0.150 

-0.243 

0.860 

0.845 

0.034 

-0.021 

3 

1.189 

1.083 

0.052 

-0.282 

0.905 

0.953 

0.089 

-0.001 

4 

1.0911 

0.957 

0.011 

-0.261 

0.951 

1.014 

0.117 

0.002 

8 

0.955 

0.911 

-0.014 

-0.114 

1.020 

1.043 

0.088 

-0.026 

OO 

1.0 

1.0 

0 

0 

1.0 

1.0 

0 

0 


Table 7. Stress intensity factors for a crack perpendi- 
cular to the external load (R=2a, c = a, Kt=1.6, 
k 2 = 1.8) . 



y 2 = ° 


y 2 = 23lJ l 


b/a 

k ll 

k 12 

k 21 

k 2 2 

k ll 

k 12 

k 21 

k 22 

/3/2 

2.8 

1.562 

. 

->• 00 
2.700 

0.206 

— OO 

-0.790 

0.781 

0 

0.483 

0.002 

-> 0 
0.167 

3.0 

1.430 

1.970 

0.177 

-0.316 

0.817 

0.613 

-0.005 

0.057 

3.5 

1.274 

1.548 

0.134 

-0.047 

0.878 

0.752 

-0.012 

-0.047 

4 

1.193 

1.364 

0.107 

0.023 

0.914 

0.833 

-0.012 

-0.068 

6 

1.070 

1.109 

0.048 

0.032 

0.970 

0.952 

-0.006 

-0.041 

8 

1.036 

1.049 

0.025 

0.015 

0.985 

0.980 

-0.002 

-0.021 

10 

1.022 

1.028 

0.014 

0.007 

0.991 

0.989 

-0.0004 

-0.012 

CO 

1.0 

1.0 

0 

0 

1.0 

1.0 

0 

0 





























Table 8. Stress intensity factors for a crack perpendi- 
cular to the external load (R = 2a , c = 2 . 2a , <-^ = 1.6, 

<2 = 1 • 8 ) . 


2 0 p 2 23y l 


b/a 

i — i 
rH 

44 

k 12 

k 21 

K 22 

k ll 

k 12 

k 21 

k 22 

0 

0.146 

0.146 

-0.113 

0.113 

1.199 

1.199 

0.127 

-0.127 

0.5 

0.549 

0.008 

-0.124 

-0.101 

1.096 

1.214 

0.202 

-0.026 

1.0 

0.952 

0.196 

-0.043 

-0.372 

0.972 

1.156 

0.212 

0.031 

2 

1.315 

0.926 

0.092 

-0.526 

0.847 

0.994 

0.122 

0.035 

4 

1.188 

1.268 

0.123 

-0.060 

0.911 

0.871 

0.011 

-0.066 

6 

1.084 

1.126 

0.0765 

0.025- 

0.961 

0.941 

-0.0016 ' 

-0.060 


Table 9. Stress intensity factors for a crack perpendi- 
cular to the external load (c=R, b-a = 0.2R; a= constant, 
R variable, or R = constant, a variable, tc^=1.6, kt 2 = 1-8). 


2 23y l 



R/a 

k ll 

k 12 

i — i 
CN 
44 

k 22 

k ll 

k 12 

k 21 

k 22 

5 

0.366 

0.058 

-0.102 

-0.021 

1.120 

1.096 

0.0256 

0.124 

3 

0.882 

0.257 

-0.0774 -0.349 

0.984 

1.063 

0.139 

0.150 

2 

1.247 

0.600 

0.0935 -0.826 

0.868 

1.022 

0.172 

0.162 

1.5 

1.354 

0.926 

0.246 

-1.226 

0.829 

0.990 

0.168 

0.162 

1.0 

1.312 

1.384 

0.400 

-1.719 

0.839 

0.951 

0.165 

0.153 

0.75 

1.225 

1.679 

0.453' 

-1.972 

0.866 

0.929 

0.176 

0.144 

0.5 

1.092 

2.180 

0.486 

-2.239 

0.899 

0.909 

0.194 

0.129 


(c) The Anti-Plane Shear Problem for 
Arbitrarily Located Crack 

Similar results for the anti-plane shear problem 
are shown in Figures 4-9 (for the geometry and notation 
see the inserts in the figures and note that in these 
figures b is the distance from the inner crack tip to 
the center of the inclusion) . k(+a) shown in Figures 

4-6 are defined by 

k(a)p Q /a = k 3 (a) = lim /2 (x-a) ^ (x,c) , 

x-*a y 

k(-a)p /a = k, (-a) = lim /2 (x+a) t, (x,c) , (18) 

x-*~ a 
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and in Figure 7 we have 

k(+a) = k 3 (+a)/(p Q /R) , (19) 

where P 0 = T lyz / Y ** +0 ° • external load. In Figure 

7 the asymptotic slopes for (a/R) -> 0 are 0.57, 1.0, and 
1.47 for m = 0, m=l, and m=23.3, respectively, where 
m=M 2 /l J l- These numbers are obtained from the solution 
of an infinite plane with a central crack subjected to 
crack surface tractions given by Figures 8 and 9 at 

the appropriate locations (which give the stresses in 
the uncracked matrix containing an elastic inclusion or 
a hole) . The details of the anti-plane problem may be 
found in [20] . 

As mentioned earlier, the general three-dimensional 
matrix- inclusion-crack problems seem to be analytically 
intractable. However, under certain simplifying assump- 
tions, the elasticity solutions for a group of three- 
dimensional problems can be obtained. As an example we 
may mention the problem of an elastic matrix containing 
a penny-shaped crack and reinforced by elastic filaments 
which are distributed perpendicular to and around the 
crack. The solution of this problem is given in [21]. 
Some sample results taken from [21] may be found in [22]. 
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